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Abstract 

In connection with the investigations of initial stages of appearance 
of turbulence in the current-carrying mediums and also the investigations 
of relaxation oscillations in thin-film bridges of high-temperature super- 
conductor Y Ba2Cu2,OT- x some problems of estimation of asymptotical 
parameters of chaotic dynamic systems are considered. The mechanism 
of measuring modes matching (discretization by time and by amplitude 
at the limited buffer memory of registering device) to the fixed procedure 
of signal handling is discussed. It is shown that the problem of estimation 
of embedding dimension can be reduced to the investigation of sufficiently 
simple phase transitions at the integer- valued finite-dimensional lattices. 

1 Introduction 

As number of experimental works devoted to noncquilibrium phase transitions 
is increasing, the problem of obtaining of asymptotical estimations of natural 
chaotic systems parameters becomes more urgent. 

As a rule, records of long time series are considered to solve this problem 
(see review [0). 

This work was accomplished in connection with Q , || where an initial stage 
of appearance of turbulence at the electrical explosion of conductors is investi- 
gated, and also in conection with || , |5) devoted to experimental investigations 
of relaxation oscillations of current in the thin film of high-temperature super- 
conductive ceramic. 

One of most important parameters of a dynamic system in the steady-state 
mode is an embedding dimension — number of essential variables of the problem 
(dimension of a half-stream). This value is of special interest for the theory too, 
in spite of that for many nonlinear problems of mathematical physics (including 
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boundary problems of magnetohydrodynamic, see review ||). As always, the 
problem of choice of number of noncquilibrium modes that should be kept at the 
transition from the base system of differential equations in partial derivatives 
to the final system of ordinary differential equations remains important. 

Almost all practical handling algorithms are based in one way or another 
on the classical result 0, according to which the necessary estimation can be 
made having the signal record from only one of stream variables. 

We also will use this result understanding that the Takens theorem is valid 
only for a smooth mappings, but in the experiment we deal with the signal 
discretizated by time and by amplitude (by the space) as well. 

It is well known (see for example review |h]]) that the small perturbations 
of unstable dynamic systems can cause qualitative change of their behaviour. 
For example, numerical modeling isn't able to answer some questions because 
of rounding errors. As to unstable cycles, the amplitude discretization can lead 
to their stabilization and the time discretization — to their destruction. 

Besides, we will actively use the construction that (as we know) is first used 
by Gregory E. Falkovich in the article || (see also ||), devoted to experimental 
estimations of embedding dimension and fractal dimension of chaotic attractors, 
appearing at the initial stage of hydrodynamic turbulence in the Couette flow. 

We will dwell on this construction in more details noting only, that there is 
used the almost obviously fact: if the next after the marked (reference) zone of 
signal single reading is functionaly connected to the previous reading ^, than 
meeting the similar zone of signal we can hope that reading immediately follow- 
ing the referenced one is similar to single one. 

In the first part of this work the description of experiment is presented and 
the methodical questions concerning to the estimation of embedding dimension 
by experimental data are discussed: accounting of finite accuracy of signal am- 
plitude measurements makes impossible the direct application of known Takens 
algorithm. 



2 Experiment 

The physical results obtained in experiment presented partly in Q , Q . We will 
discuss only methodological side of it. The relaxation oscillations of current and 
voltage in the thin film high-temperature superconductor bridges Y B<i2Cu?,0-j~x 
(critical temperatures — T c = 86 — 88 K, density of critical current — j c = 
10 5 — 10 6 A/cm 2 at 77 K and absence of external magnetic field) were inves- 
tigated. The samples were fabricated by direct current magnetron sputtering 
with target of stehiometric consistant onto the substrates made of monocrystals 

3 In |Q, |?j the three-mode approximation is used and the theoretical and the experimental 
arguments in favour of this low-mode model are presented. 

4 It take place if the number of readings in the reference zone is equal or greater than the 
number of independent variables of dynamical system 
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SrTiOz ZrOziy-iOz). The dimensions of bridges were: thickness — ~ 0.3 /xm, 
width — 0.7 - 1.0 mm, length — 1-7 mm. 



Figure 1: Schematic diagram of measurements. At the obtaining of steady-state 
voltage-current characteristic Lq = 0. 

The sample was connected to the electrical circuit according to the schematic 
diagram shown in the Fig. 1 and was plunged into the dewar with the liquid 
nitrogen. The measurements of current and voltage was made by using standard 
four-points scheme. The oscilloscope C9-8 was used for registration of oscillo- 
gramms. It has 8-digits analog-digital converters and the buffer size for signals 
recording is 2048 bytes. The steady-state voltage-current characteristic shown 
in Fig. 2 were obtained in the mode with fixed source voltage and the connect- 
ing the sample into the circuit in series with load resistor (without inductance) . 
The portion of voltage-current characterictic wit negative differential resistance 
shown in the Fig. 2 defines development of instability ||. The oscillogram of 
voltage in circuit with inductance (L ^ 0) is shown in Fig. 3. It is obvious 
that we deal with well expressed periodic mode. 

Note that current I(t) is directly among of the set of phase variables govern- 
ing the system dynamic. It comes in the form of the integral of j p — laminar 
consistant of current taken over cross-section of superconductive film (see Q, 



The voltage U(t) — J Eidx, where Ei is the longitudinal component of elec- 



trical field on the film's surface, I is the distance between the measuring elec- 
trodes. 

In the above considerations I(t), U(t) are continuous and rather smooth 
functions of time. 

The analog signals I(t), U(t) are applied to the input of measuring device 
that makes encoding (discretization by amplitude) of both signals with the 
regular time step r. On completing of one measuring cycle we have the ordered 
array of code pairs, filling in the buffer memory of oscilloscope. Every code can 
be interpreted as integer positive number in the range [0, 2 m — 1], m = 8. It is 
possible to assume that this number represents the result of integer division of 



Figure 2: a) The steady-state voltage-current characterictic; b) The dynamic 
voltage-current characteristic in periodic mode (Lq ^ 0). 
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Figure 3: The oscillogram of voltage in HTSC film. The selected fragment 
shown in Fig. 6. 



Figure 4: The normalized dynamic voltage-current characteristic in the periodic 
mode. The selected fragment is shown in Fig. 5. 

"instant" value of corresponding measuring signal by the interval between the 
levels of it's discretization. 

Almost all possible range of codes is usually in effect. At these conditions 
the aprioral characteristic of accuracy of measurement of signals by amplitude 
(by space) is the value e = 2~ m Q At e — ^ we can assume that the process of 
measurement reduces to time discretization of continuous signals. 

The so called dynamic voltage-current characteristic of film is shown in Fig. 
4. Every it's point corresponds to the pair of codes (I, U) obtained at the same 
moment. This figure can be considered as a point mapping of a segment of 
phase trajectory of the system onto the plane of measurements. The large- 
scaled fragment of Fig. 4 is shown in the Fig. 5. We can see that the points of 
trajectory of discretizated system are located at the nodes of the integer-number 
mesh (lattice). 

One more effect of space discretization is the appearing of the consecutive 
identical codes in sequences of codes (flat portions in oscillograms) (see Fig. 6). 

What hinders the situation is that the perturbations effected by the mea- 
surement procedure are not lone in the system. 

On amplitude measurements. The amplitude of signal and the average point 
of code scale (the value of code at the zero signal) we select so that the realization 
will not contain codes located at the border of the scale. On the choice made the 
system trajectory calls upon only the internal points of coding space. If to > 8 
this choice does not lead to significant loss of accuracy. The corresponding 
component of error is as a rule in the range 



On choice r — time step value. The limitation of buffer memory size essen- 
tially influences this choice. 

Too fine quantization will cause the uneffcctive filling of N volume by the 
long sequences of identical repeating codes in limiting case by single code. To 

5 The situation like this (so called regular e — discretization) appears due to the errors of 
rounding at the numerical modeling of dynamic systems JlOl 




Figure 5: The fragment of Fig. 4. 
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Figure 6: Fragment of the oscillogram of voltage from Fig. 3. 



Figure 7: The model signal d = 3, n = 8. 



make r greater than characteric correlative time is uneffective too. 

In actual practice, at least at the initial stage of investigation the principle 
of "peer interest to space and to time" is in effect and the step magnitude To 
is set to make the measurement error of characteric time interval (in this case 
— of period T) to be close to e. For example, if e = 2~ 8 , N = 2 11 bytes than 
on setting r in accordance with aforementioned principle it will be registered 
K ~ Ne, Ko ~ 8 periods of oscillations (See Fig. 3). 

To check if it is possible to consider the trajectory as a cycle it is necessary 
to have K min > 2, i. e. r m i„ > 0, 25ro- One can evaluate K max specifying the 
limit value of error in the estimation of period length. So, for e max = 3%, we 
have K max < 32, i. e. < 4r . 

3 Handling of data 

At the handling of experimental data we will try not to go off the frames of 
integer numbers set that is inherent to this data. In some cases it makes pos- 
sible to reduce the problem to investigation of rather simple structural phase 
transitions in the integer-number finite-dimensional lattices. The benefits of 
this method, we think, make up for possible inconveniences. In connection to 
this circumstance we will consider the main procedures in some simple model 
examples. 

3.1 Integer- number model 

Let us assume the model signal fragment of which is shown in Fig. 7: 



The first step envisages the building from (1) vectors (corteges) with the 
length k by the next algorithm: 
i) Fixing initial vector 



U(t) = {01212321012123...} 



(1) 



U(t) = {Uj}, j 1 11 V. 




<j+i 



(2) 



5 



j = l,...,N-k. 

As a result, the sequence appears with length of (N — k) fc-components 
vectors. In our example for k = 1 it is just starting sequence, for k = 2 — 13 
pairs: 

< (01)(12)(21)(12)(23)(32)(21)(10)(01)(12)(21)(12)(23) >, (4) 
For k = 3 — 12 triplets. 

< (012)(121)(212)(123)(232)(321)(210)(101)(012)(121)(212)(123) > (5) 

It is necessary to find k — d so that the cortege < £f , > is a cycle of 

period n with the following properties: 

a) there are no identical vectors in the limits of a period 

tf^, i^J, *,i<n; (6) 

b) the vectors separated one from another with period exactly are identical 
ones 

€f+nK = (7) 

here K is the ordinal number of a period. 

As it was specified, the sequence (4) is not a cycle because it contains re- 
peating pairs. Nevertheless, (5) is a cycle of period n = 8. 

The desired embedding dimension is U(t) d = 3. Note, that here d is integer- 
valued. It should be considered as a right border of noninteger- valued dimen- 
sions that can characterize the dynamics of system in the chaotic modes of its 
functioning. 

It is clear that in our case one has to understand the identity (nonidentity) 
in a certain sense. Before to do it, let's recall that (i + k + l)th reading of main 
sequence Ui + u+i directly follows £ 4 fc , and Uj + k+i follows £j\ 

Let us settle some criterion of vectors dist(£f,£j) resemblance and 

assume the small e, generally speaking, depending on i, j. The criterion that on 
k > d (k + l)-th reading is a function of k previous readings can be formalized 
in the following way. 

If on k > d there exists an e- vicinity of zero (e p , e r ) so that from 

dist(tf,tij)<e P (i,j) (8) 
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Figure 8: Meshes for the signal in Fig. 7; k = 1, I = 1 -j- 4. 
Figure 9: Meshes for the signal in Fig. 7; k = 1, Z = 5 -r- 8. 

follows 

dzst(t/ i+ fe + i,U J - +fe+ i) < Cr(i,j), (9) 
then fc = d is an desired evaluation. 

In cited above || Euclidean distances are set as a degree of proximity 

p>< = dist(l;)=\\$-$\\, (10) 

and 

r fc = dist{U) =|| C/ i+fe+1 - U j+k+1 ||, (11) 

where 1 1 • 1 1 is the norm of corresponding differences. It is not very convenient 
because p k , r k appear to be real numbers albeit initial data belongs to the field 
of integer numbers. According to the assertion at the beginning of section 3, 
we will use a module (not a norm) of corresponding differences as a degree of 
proximity (unless specified otherwise): 

p k =\tf-$\, (12) 

and 

r k =| U i+k+ i-U J+k+1 | . (13) 

Using the algorithm (1), the algorithm i)-ii), and also (12)-(13) we build the 
plots r k (p k ) consecutively comparing the selected pair < ^ k ^Ui+k+i > to the 
rest of pairs j > I. 

This procedure makes a mapping of the initial trajectory (1) at k = 1 into 
the square | 3 X 3 | , at k = 2 — into rectangle | 3 X 6 | , at k = 3 — into | 3 x 9 | , 
and so on. At the first step j — I, so that the trajectory (p, r) starts from the 
point (0,0) and on completing the cycle comes back to this point. 

Varying I — 1, n for each k = 1, 2, 3 we get n meshes showed in Fig. 8 - 
17. From the point of view of topology all meshes represent multy-beam stars 
with centre in the point with coordinate (fc, 1). 

Let choose as a 



Figure 10: Meshes for the signal in Fig. 7;fe = 2,i = l-=-2. 



7 



Figure 11: Meshes for the signal in Fig. 7; k = 2, I = 3 ^ 4. 



Figure 12: Meshes for the signal in Fig. 7; k = 2 



I = 5 -=-6. 



e p = min{p ,ps} 



(14) 



where po = ro is abscissa of the centre of the star at k = 1, and ps evaluate 
as a point of intersection of the line r = p, drawn from the point of contact 
of the plot with the axis of ordinatcs to the point of intersection with its top 
delimiting border (see Fig. 8). For this model ps = ro = 1. 

Beginning from the Fig. 8 - 9 (k = 1) in two cases (I = 1,6) of eight the 
structure of mesh meet the demands of criterion (8), (9). Namely, in the area 
e p < 1 all its points do not fall beyond the border line r = p. 

In Fig. 10 - 13 (k = 2) four meshes from eight possess this property, in 
Fig. 14 - 17 (k = 3) all of eight. As to the rest of meshes in Fig. 8 - 9, 10 
- 13, presence of the "bad" node with coordinate (0,2) within their structure 
prevents the meeting of criterion (8), (9). 

Thus the structural phase transition in a discrete finite- dimensional lattice 
is present. The subset of meshes (number of meshes n a ) whose structure com- 
prises the nodes located beyond the line r = p acts as initial stage, np meshes 
represents the new /3-stage and in the sum n a + np = n where n is a period of a 
cycle. For this phase transition k act as the parameter of order and k — d is its 
critical value. Note two more peculiarities: a) the distance between the location 
of "bad" node (0, 2) in the initial stage and "good" node (2, 2) in the new stage 
is constant along the lattice (gap) and equal 2; b) node with coordinate (0, 1) 
is always empty. 

The model signal is an ideal in many respects. For example, if we take (1) 
and trace the dynamics of volume of elementary cell of phase space 



I = 1, n, ...TV — 3, we will come to the conclusion that AV(l) is a unit cube 
propagating without deformation. 

There is no additive noise in (1) and the amplitude transitions occur only 
between two contiguous levels. As a result in all of plots the centre of the star 
is located at the ordinate closest to the zero ro = 1. In general, r > 1. 
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(15) 



Figure 13: Meshes for the signal in Fig. 7;fc = 2,Z = 7-^8. 
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Figure 14: Meshes for signal in Fig. 7; k = 3, I — 1 -f- 2. 

Figure 15: Meshes for signal in Fig. 7; fc = 3, I = 3 -f- 4. 

There never appear two (and more) identical sequential codes in the model 
example (1) but it is by no means nearly so in practice. For example, Fig. 3 
contains (see Fig. 6) long (up to 32 elements) sequences of identical codes — this 
circumstance can radically deteriorate the quantity of estimation d according to 
this algorithm. 

As illustration to the last assertion we consider the sequence: 

{011221122332211011...}, (16) 

obtained from (1) by repeating twice of all codes besides 0. The meshes for 
(16) are shown in Fig. 18. The following changes occured as against Fig. 8 - 
17: 

— it is no longer a star with well emphasized centre; 

— the loops corresponding to the consecutive identical pairs (p, r) appeared 
into the graph structure besides simple cycles and short sequences; 

- the trajectory repeatedly visits the axis of ordinates at the point (0, 1) 
that did not occure before. Nevertheless, the phase transition takes place at 
k = 3. 

The sitution changes drastically when the maximum length of a sequence of 
identical codes in signal Al max > d. From this moment the Al max governs the 
estimation of d. Fig. 19 presents the mesh (k = 3, 1 = 2) for the signal that 
differs from (16) by adding of third unity into the sequence of unities following 
the zero. The mesh in Fig. 19 is shifted by one step to the left as against 
Fig. 15. This has the effect of shifting the good node (2,2) to the point with 
coordinates (1,2) converting it to the bad one. 

Note some effects deteriorating the quantity of estimation of d due to in- 
creasing systematic part of error. With the increasing of k the space taken by 
the plot is also increasing (number of nodes of the coordinate mesh grows) . On 
fixing: (1) volume of realization and (2) size of zone close to zero p = [0, e p ] that 
we are interested in, the density of depicting points decreases. 

As illustration we consider this effect on example of sequence N = 2048 
bytes generated by the pseudo-random number generator for k = 1,3,6 (Fig. 
20). While at k = 1 depicting points evenly cover the space of the plot, on 
increasing of k they contract into the narrow band in its centre. []. 

6 One can "lose" rarefied zone near the zero starting calculation r(p) from j = I (not from 



Figure 16: Meshes for signal in Fig. 7; k = 3, I — 5 -j- 6. 



Figure 17: Meshes for signal in Fig. 7; k = 3, / = 7 -f- 8. 

Figure 18: Meshes for sequence (16). 

Completing the subsection (3.1) we address again to J|] (see Introduction) 
where the concept of uppermost envelope for the plots (p, r) is adopted and its 
monotonicity near the zero is stated. Actually, in the limit e — > one may use 
term envelope r(p). At k > d in the zone of a small p < e p scales from 

/=IUf-^IHo, (17) 

follows 

r d =|| U{t l+d+1 ) - U{t j+d+1 ) ||-> 0. (18) 

In real life e ^ 0; the sequences of repeating codes may be present in sig- 
nal and the longest of them directs eventually the quantity of estimation of 
embedding dimension by Takens-Falkovich algorithm. 

4 Conclusions 

Thus, we showed that signal amplitude discretization makes it impossible the 
direct application of classical Takens-Falkovich algorithm for the direct estima- 
tion of embedding dimension. One can note that various signal sections are 
disparate in the sense that they are able to yield the different estimations for 
d. This makes us to choose the fixed (depending on signal properties, accuracy 
of its amplitude and realization length measuring) time discretization step and 
to solve the phase transition problem onto the lattice (in accordance with the 
algorithm presented in the section 3). 

In the second part of this paper we will consider some techniques of obtaining 
the desired estimation by the example of Lorenz system (operating in various 
modes) and we are going to present the results of processing of experiment 
described in section 2 obtained by means of program taught by Lorenz system. 



j = 1 + 1). In this case (using standard procedures of of plot drawing), the absciss of the left 
low corner will be defined by minimum p ^ in the particular realization. 



Figure 19: Meshes for signal in which as opposed to (16) the third unity is added 
to the sequence of unities following the first zero. 
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Figure 20: The illustration of systematic error stipulated by finiteness of real- 
ization length N. 
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